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ABSTRACT 

Context. The spectrum of cosmic rays (CRs) is affected by their escape from an acceleration site. This may have been 
observed not only in the gamma-ray spectrum of young supernova remnants (SNRs) such as RX J1713.7— 3946, but 
also in the spectrum of CRs showering on the Earth. 

Aims. The escape-limited model of cosmic-ray acceleration is studied in general. We discuss the spectrum of CRs running 
away from the acceleration site. The model may also constrain the spectral index at the acceleration site and the ansatz 
with respect to the unknown injection process into the particle acceleration. 

Methods. Analytical derivations. We apply our model to CR acceleration in SNRs and in active galactic nuclei (AGN), 
which are plausible candidates of Galactic and extragalactic CRs, respectively. In particular, for young SNRs, we take 
account of the shock evolution with cooling of escaping CRs in the Sedov phase. 

Results. The spectrum of escaping CRs generally depends on the physical quantities at the acceleration site, such as 
the spectral index, the evolution of the maximum energy of CRs and the evolution of the normalization factor of 
the spectrum. It is found that the spectrum of run-away particles can be both softer and harder than that of the 
acceleration site. 

Conclusions. The model could explain spectral indices of both Galactic and extragalactic CRs produced by SNRs and 
AGNs, respectively, suggesting the unified picture of CR acceleration. 

Key words. Acceleration of particles - ISM: cosmic rays - ISM: supernova remnants - Galaxies: jets 



1. Introduction 

The origin of cosmic rays (CRs) has been one of the long- 
standing problems. The number spectrum of nuclear CRs 
observed at the Earth, M(E) cx E~ s , shows a break at 
| the "knee" energy (w 10 15 5 eV), below which the spec- 
tral index is about s w —2.7 (|Croninlll999f ). Because of the 
energy-dependent propagation of CRs, the spectral shape 
at the source is different from that observed at the Earth. 
Taking into account the propagation effect, the source spec- 
tral index has bee n well constrained as s _ ~ 2.2 — 2.4 in var- 
ious models (e.g.. IStrong fe Moskalenkolll998t [Putze et alj 
2009). This value of s has been also inferred in order 
to explain the Ga lactic diffuse gamma-ray emission (e.g., 
IStrong et al] |2000). This fact may give us valuable insights 
on the acceleration mechanism of CRs. 

Mechanisms of CR acceleration have also been studied 
for a long time and the most plausible process is a diffusive 
shock accelerati o n (PSA) (iKrvmskvl Il977t lAxford et aTl 
fl977t iBeil U97l iBlandford fc Ostrikerl [l978h . Very high- 



ber o f nuclear particles are accelerated (jDrurv fc Volkl 
119811 : iMalkov fc Drurvl l200lh . There are several observa- 
tional facts which are consistent with the predictions of 
such nonline ar model (jVink fc Laming 2003; Bamba et alJ 



20031 l2005aira: iWarren et all 120051 : lUchivama et all [20071 : 
120091) . predicts, however, the 



Helder et al 



energy gamma-ray observations have revealed that the ex- 
istence of high-energy particles at the shock of young su- 
pernova remnants (SNRs), which supports the DSA mech- 
anism as well as the paradigm th at the Galactic CRs 
are produced by young SNRs fe.g.. [En omoto ct "aTl [20021 : 
lAharonian et alj 12004 120051 iKatagiri et all 120051 ). Recent 
progress of the theory of DSA has revealed that the back- 
reactions of accelerated CRs are important if a large num- 



harder spectrum of accelerated particles at the shock than 
f{p) cx p~ 4 (corresponding to s = 2) where p is the momen- 
tum of CRs, in particular, near the knee energ y (jMalkovl 
[l997tlBerezhko fc Ellisonll 19991: iKang eTldlf^OOlf) . This fact 
apparently contradicts with the source spectral index of 
s w 2.2 — 2.4 inferred from the CR spectrum at the Earth. 
Even in the test-particle limit of DSA, such a soft source 
spectrum requires a shock with the small Mach number 
(M 10), which is unexpected for young SNRs. 

There are several models of DSA, depending on the 
boundary conditions imposed. Different models predict dif- 
ferent spectra of CRs dispersed from the shock region. So 
far, the age-limited acceleration has been frequently con- 
sidered as a representative case (§ [2]). In this model, all 
the particles are stored around the shock while accelerated. 
When the confinement becomes inefficient, all the parti- 
cles run away from the region at a time. Then, the source 
spectrum of CRs which has just escaped from the accelera- 
tion region is expected to be the same as that at the shock 
front. Therefore, this model predicts that the source spec- 
trum is the same as that of accelerated particles, which is 
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typically harder than the observed one. In this paper, we 
consider an alternative model, the escape-limited acceler- 
ation, to explain the observed CR spectrum at the Earth 
C § EJ) - This model is preferable to the age- limited accelera- 
tion when we consider observational results for young SNR 
RX J1713.7— 3946, of which TeV 7-ray emission is more 
precisely measured than any others (§ 14. ip . 

The nature of CRs with energies much higher than the 
knee energy is also still uncertain. While CRs below the 
second knee (~ 10 175 eV) may be Galactic origin, the 
highest energy CRs above ~ 10 18 ' 5 eV are believed to be 
extragalactic. Po ssible candidates are active galactic nu- 
clei (AGNs) (e.g-jBiermann fc Strittmatterlll987t[Takaharal 
ll990HRachen fc Biermann|1993HPe'er et alj200jj 
ray bursts (IWaxmanll 19951: IVietrilll995| 
magnetar s (lArons 2003; M urase et al.ll2009l 
galaxies (|Kang et all Il997t llnoue et all 1200 
mediate energy range from ~ 10 17 - 5 eV to 
more uncertain 



cles, pm S °' , is determined by i ac 
are the age of the shock and the acceleration time sc ale, re- 
spect ively. When we consider DSA, t acc is given by (jPrurvl 
1983) 



i ag0 , where t age and t acc 



amma - 
Muraseetaljl2006), 



3) an d clusters of 
007t ). The inter- 
lO 18 ' 5 eV is 

Both the Galactic and extragalactic ori- 
gins are possible and it may just a trans ition between the 
two. As the extragalactic origin. AGNs jBerezinskv et al.l 
l2006t ) . clusters of galaxie s (jMurase et al.ll2008l ) and hyper- 
novae (|Wang et al.ll2007l) have been proposed so far. 

Among these possibilities, AGN is one of the most 
plausible candidates for accelerators of high-energy CRs, 
because it can explain the UHECR spectrum above ~ 
10 17 ' 5 eV assuming the proton composition. In such a pro- 
ton dip model, the source spectrum of UHECRs is Af oc E~ s 
with s = 2.4-2.7, depending o n models of the source evo- 
lution (jBcrczinsk v et all 120061 ). The required source spec- 
tral index of s = 2.4-2.7 can be explained by several pos- 
sibilities. First, it can be attributed to the acceleration 
mechanism it self. One can consider non- Fermi acceleration 
mechanisms ([Berezinskv et all 12006) or the t wo-step diffu- 
sive s hock acceleration in two different shocks (jAloisio et alj 
120071 ). Second, the index can be attributed to a superposi- 
tion of many AGNs with different maximum energies, and 
one can suppose that AGNs with di fferent luminosities may 
have different maximum energie s ([Kachelriess fe Semikoa 
I2006T) . Recently, iBerezhkol (|2008D proposed another possi- 
bility under the cocoon shock model. In this cocoon shock 
scenario, different maximum energies can be interpreted as 
maximum energies of escaping particles at different ages of 
AGN jets. Although it is very uncertain whether efficient 
CR acceleration occurs there, this scenario would also be 
one of the possibilities to be investigated in detail. 

The organization of the paper is as follows. After the 
brief introduction of the age-limited model of the CR ac- 
celeration (§0, we study the escape-limited model in §[3] 
For a simple understanding, the general argument in a sta- 
tionary, test-particle approximation is given in § 13.11 Then, 
we derive the formulae of the maximum energy of accel- 
erated particles in § 13.21 and of the spectrum of escaping 
particles in § 13.31 We consider the applications to young 
SNR and AGN in § g] and § respectively. Section [5] is 
devoted to a discussion. 



2. Maximum Attainable Momentum in the 
Age-limited Acceleration 

For comparison with the escape-limited acceleration, we 
briefly summarize the case of the age-limited acceleration. 
In this case, the maximum momentum of accelerated parti- 



al - U 2 



Di( P ) , D 2 (p) 



u 2 



(1) 



where D(p) and u are the diffusion coefficient as a function 
of the momentum of accelerated particles and the velocity 
of the background fluid, respectively. Subscripts 1 and 2 
represent upstream and downstream regions, respectively. 
For simplicity, we assume the Bohm-type diffusion, i.e., 



D 1 {p) = D 2 (p) 



3eB 



(2) 



where B, r) g and m p are the magnetic field strength, the 
gyro-factor and the proton mass, respectively. Taking into 
account that the fluid velocities are related to the shock 
velocity, u s h, as u% = au 2 = u s h, wh ere a is the compres- 
sion r atio of the shock, we derive (e.g.. lAharonian fe Atovanl 
119990 



n( a S C ) 



1 eB 



a(a + 1) 



sh'-agc 



(3) 



3. Escape-limited Acceleration 

In the frame work of DSA, accelerated particles are scat- 
tered by the turbulent magnetic field, and go back and forth 
across the shock front. Upstream turbulen ce may be e xcited 
by the accelerated particles themselves (jBelll Il978t ). and 
the magnetic field strength of s uch turbulence is t heoreti- 
cally expected to be strong fe.g jLucek fe Beil2000( ). There 
are observational evidences suggesting that CRs are re- 
sponsible for substantial amplification of the ambient mag- 
netic field in the precursors of shock fronts in SNRs, and 
that such magnetic tu rbulence well confines the particles 
around the sh o ck front (|Vink fc Lam ing 2003: lBamba et all 



20031. l2005allbl: 



Uchivama et al 



Yamazaki et all 12004 iParizot et all 12006 



20071 ). leading to the efficient CR accel- 



eration. 

The spectrum of accelerated particles is affected by 
the spatial and spectral structures of the magnetic tur- 
bulence through the process in which the particles es- 
cape from the shock toward far upstream region. There 
are mainly two scenarios of the escape model considered 
so far; one causes the effect on the boundary in the mo- 
mentum space, and the other causes the effect on the spa- 
tial boundary. The former comes from significant decay of 
the wave amplitude below t he wave number fc m j n of the 
spect rum of the turbulence ([Reynolds! 119981 : iDrurv et alj 
2009) . Particles with the Lorentz factor above 7 C satisfying 
the resonance condition, u — k mm v — ^c(7c) = 0, where 
fi c (7) = eB/^m-pC is the cyclotron frequency, are not con- 
fined around the shock front and escape into far upstream 
region. In t his context, the escape flux was calculated prev i- 
ously (e.g. JPtuskin Zirakashviiil2005l : lDrurv et al.ll2009t) . 
The latt er effect has been recently discussed by several 
authors (IPtuskin fc Zirakashvilil [2001 iReville et all 120091 : 
ICaprioli et al.l [2009f ). The turbulence generation may be 
connected with the flux of accelerated particles themselves. 
Hence, in the region far from the shock front, the flux of 
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high-energy particles is small and wave excitation is less sig- 
nificant. If the accelerated particles reach the region, they 
are dispersed into the far upstream region. Let £ be the dis- 
tance from the shock beyond which the amplitude of the up- 
stream turbulence becomes negligible. Characteristic spa- 
tial length of particles penetrating into the upstream region 
is given by D(p) /u s h- As long as D(p) /u s h <C £, the particles 
are confined without the significant escape loss, and they 
are accelerated to higher energies. On the other hand, when 
their momentum increases up to sufficiently high energies 
satisfying D(p)/u s h *S £, their acceleration ceases and they 
escape into the far upstream. Therefore, the maximum mo- 
mentum of accelerated particles in this scenario is given by 
the condition D(p)/u s h ~ £■ In the following of the paper, 
we consider the escape-limited model, where the maximum 
energy is essentially determined by D(p)/u s h £■ 

3.1. A simple case of stationary, test-particle approximation 

In order to take an essential feature of the escape-limited 
acceleration, we calculate the escape flux and the maxi- 
mum attainable e nergy of accelerated particles in the sim- 
plest case (see also lCaprioli et al1l2009D . Let us consider the 
stationary transport equation 



u{x) 



dx 



d_ 

dx 



D(p) 



df 
dx 



p du df 
3 da; dp 



(4) 



with the boundary condition, f(x — —£) = 0, where x = — I 
(£ > 0) is the upstream escape boundary. The fluid velocity 
u{x) is given by 



u{x) 



ui (x < 0) 
u 2 (x > 0) 



(5) 



where u\ and ui are constants. The solution to the trans- 
por t equation in the tes t-particle approximation is derived 
as (|Caprioli et al.ll2009t) 



f(x,p) = f (p)- 

where fo{p) = f(x = 0,p) is given by 

dlogp' 



exp 


_D(p)_ 




exp 




ml 
D(p)_ 


1 — exp 


ml 
D(p)_ 





(6) 



fo(p) = if exp } -q 



ml 



D(p') 



'Pinj 1 — exp 

and q = 2>u\/{u\ — u 2 ). The escape flux at x = — £ is 

df 



(7) 



Hp) 



«1 f\x=-l~ 

utfoip) 



D(p) 



Ox 



1 — exp 



D(p) 



(8) 



Let us introduce a new variable y = hip and a new 
function $(?/) = ln0(p). Then, we expand <&(?/) around its 
maximum value at y = y m . To do this, we calculate the 
first and the second derivatives as 



d$ 

d V 1 - exp 



1 



mi 
D{y) 



u x £ dln£> 



D(y) dy 



(9) 



d 2 $ 
dy 2 



ml 



+ 



1 — exp 



1 — exp 



ml 



u x £ dlnD 
D(y) dy 



dln£> 

dy 



D{y) 



dln£> 
dy 



d 2 \nD 

dy 2 



.(10) 



In the following, we consider the case of Bohm diffusion, 
D{p) = Dq(p I 'rripc) . Then, one can find 



Pn 



ui£m p c 
qD 



(11) 



(where d$/dy — 0) so that we obtain 

1<9 2 $ 
2~dy- 



Hv) = *(2/n 



, 2 "(ym)(y - 



+ 



= In 



Ml/o(Pm) 



1 - e~i 



(12) 



where ^ 2 = 2(1 — e~ q )/q < 1. Note that the quantity p m 
given by Eq. (|lip also plays the role of maximum momentum 
of the accelerated particles at the shock (x = 0) because one 
can see f Q (p ) oc p~ q for p < p m while f (p) cx exp(-p/p m ) 
for p m < p ()Caprioli et al.ll2009l ). 

Finally, going back to the function </>(p), we obtain 



UlfoiPm) 

1 - e~i 



exp 



lnp — In p n 
I 



+ 



(13) 



One can clearly see, from Eq. (flU)) . that particles with mo- 
mentum around p m escape the shock region most efficiently. 

3.2. The maximum energy of accelerated particles 



We have seen in § 13.11 that in the stationary, test-particle 
case, the quantity p m given by Eq. (|11|) plays the role of 
maximum momentum of the accelerated particles at the 
shock. Taking this fact into our mind, we assume that in 
the more general escape-limited case, the maximum mo- 
mentum, pm SC \ is determined by 

„^ SC) ) , 



Ul 

Given that 



D{p) 



r) g m p c° 



3eB V m v c 



(14) 



(15) 



which is the same as in the age-limited case [Eq. ([2]) in § [2] , 
we obtain 



qr) g c 2 



(16) 



Pi 



Since q = 3a /(a — 1), we find from Eqs. ^ and (fll)]) 

(esc) D 



(ago) 



(17) 



P 

Hence, as long as £ 3> Ush^agc which is expected in the 
early phase of the shock, we don't need the escape-limited 
acceleration. However, if the shock evolves so that £ 
Ushtage, we have to consider the effect of the particle escape, 
and we cannot simply apply the well-known result of age- 
limited acceleration, Eq. ©. 
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3.3. The spectrum of CRs dispersed from an accelerator 

In this subsection, we derive the time-integrated spec- 
trum of CR particles which is dispersed from an accel- 
erator. The derivation is ess entially identical to that of 
iPtuskin & Zirakashvilil (j2005f ). However, our argument is 
simpler and more general, so the final form of the spectrum 
(Eqs. (|27f and (|28| ) is more general. Note that our for- 
malism is applicable not only to DSA but also to arbitrary 
acceleration processes. 

The proton production rate, N(jp, x)d%dp, at a certain 
epoch labeled by a parameter x, is defined as the number 
of protons with momentum between p and p + dp which 
is produced in the interval between \ and x + dx- Here x 
is the parameter which describes the dynamical evolution 
of the accelerator — it can be either simply the age, or 
the position of the shock front. It is expected that N(p,x) 
contains the term of exponential cutoff at the momentum 
Pm(x) which depends on x [see, for example, Eq. (|25[) ]. The 
number of protons with momentum between p and p + dp 
that is escaping from the accelerator at the epoch between 
X and dx is denoted by <f>(j>, x)dxdp, and we assume 



47rp 2 </>(p, x)dxdp oc N(p m (x), x)G(p, x)dxdp 



where 



G(j>, x) = exp 



lnp- lnp m (x) 



(18) 



(19) 



and £ < 1. This is because we expect that particles with 
momentum around p m (x) is most efficiently escaping from 
the source. Indeed, as shown in § 13.11 Eqs. (|18|) and (fl9|) 
are good approximation in the test-particle, stationary case. 
Here, we simply expect these assumptions are also correct 
in general. 

The time-integrated spectrum of protons which have es- 
caped at the source, N esc (p), is obtained by 



N esc (p) =4ttp 2 / <Mp,x)dx 



(20) 



In order to derive a simple analytical form, we approximate 
Eq. dHJ) as 

(21) 



G(p, x) ~ \/< • S{lnp - lnp m (x)) 



If we use a general mathematical formula for <5-functions 
for an arbitrary function g(x)' 



%(x)) 



S(x-Xo) 



(22) 



where g(xo) — 0, then Eq. (f2Tj) is rewritten as 

- s ix-p m Hp)) , 



G(p, x) oc 



p 


m(x) 






L d x _ 


X=Pm 1 (p) 



(23) 



where p m (p) is the inverse function of p m (x)i that is, 
Pm(Pm(p)) = P and P^iPmix)) = X- Using Eqs. (HU) and 
(|2"3"|) , we calculate N esc (p) as 



N esc (p) oc 



pN(p,p^(p)) 


PnMp) 


[ d X _ 


X=Pm 1 (P) 



(24) 



This is the most general analytical formula of the spectrum 
of protons dispersed from an accelerator. 

In the following of the paper, the form of N(p, x) is 
assumed to have 



N{p,x)dxdp = K(x) 



P 



x exp 



Pm(x) 



dlogxdp 



(25) 



which is a power-law with the index s and the exponential 
cut-off at p m (x)- Substituting Eq. ([25]) into Eq. (|24|) . we 
obtain 



N esc (p) oc 



P 



s K{p^{p)) 



t{p) 



<!//; 



(26) 



1\p) 



In particular, if K (x) and p m ix) are written by the power- 
law forms, such as K(x) oc x' 3 and p m ix) °c X~ Q j then, 
p^ip) P _1 ^ Q , so that Eq. (|26l) becomes 



N csc (p) oc p ' 
where 

^csc S i 

a 



(27) 



(28) 



This is the simplest form of the spectrum of CRs which are 
dispersed from the acceleration region. 

Generally speaking, in order to obtain the energy spec- 
trum of accelerated particles, time-dependent kinetic equa- 
tion should be solved. Instead, we have assumed that at 
an arbitrary epoch, the spectral form is given by Eq. (I25p . 
This assumption is justified if the spectrum at the given 
epoch is dominated by those which are being accelerated 
at that time, in other words, if the particle spectrum does 
not so much depend on the past acceleration history. For 
example, in the case of the spherical expansion, accelerated 
particles suffer adiabatic expansion after they are trans- 
ported downstream o f the shock and lose their energy (e.g., 
Yamazaki et al. 2006) , so that the contribution of the pre- 
viously accelerated particles is negligible. Strictly speaking, 
even if we consider the energy loss via adiabatic expansion, 
the energy spectrum of accelerated particles does depend 
on the past acceleration history in some cases. When we 
use the shock radius i? s h as x> the final form of Eqs. (j27|) 
and (|2"8")) is correct as long as (3 > (s — 1) (ct _1 — l) (see 
Appendix) , which is satisfied in the cases considered in § U 
and § [5j Otherwise, the form of Eq. ([23| is no longer a good 
approximation, and the final form of s csc is different from 
Eq. (J2H1) (see Appendix). 



4. Application to Young Supernova Remnants 

4.1. Inconsistency of age-limited acceleration with observed 
results of RX Jl 713. 7-3946 

RX J1713.7— 3946 is a representative SNR from which 
bright TeV 7-rays have been detected. The H.E.S.S. ex- 
periment measured the TeV spectrum and claimed that 
i ts shape was better e xplain ed by the hadronic model 
(Ah aroman et al.l I2006L I2007D . Furthermore, evidences of 
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amplified magnetic field (B > 0.1 mG) are derived from 
the width of s ynchrotron X-ray filaments (Parizot et al . 
2006; see also IVink fe Laming] 120031: IBamba et all l2003l 
2005 aH) and from time varia tion of synchrotron X-ray hot 
spots ( Uchivama et al.|[2007t^ . These facts also support the 
hadronic origin of TeV 7- rays, because the leptonic, one- 
zone emission model (e.g.. lAharonian fc Atovan1 IT999) can- 
not explain the TeV-to-X-ray flux ratio. Hence it is natural 
to assume that the TeV 7-rays are produced by the hadronic 
process althou gh there are several arguments against this 
interp retation flButtl I2008t iKatz fc Waxmanl l2008t iPlaeal 
120081) . 

In the age-limited case, Eq. ([3| reads 



cp 



(age) 



8x 



lo3 ^(3o5o^) 2 fe) TeV ' (29) 



where we adopt a = 4, and B m Q is the magnetic field 
strength in units of mG. 

H.E.S.S. observation has revealed th at the cutoff en- 
ergy of Te V 7-ray spectrum is low (jAharonian et all 
120061 I2007t ). so that in the one-zone hadronic scenario 

the maximum energy of p rotons, -Emax.p = cpm ge ^ is 
estimated as 30-100 TeV (jVillante fc Vissanil 120071 ). If 
E max , p < 100 TeV and B w 1 mG, then Eq. ([29]) 
tells us ri g > 80, implying far from the "Bohm 
limit" (r]g m 1) which is in f erred from the X-ray ob- 
servation (jParizot et al.l l2006t lYamazaki et al. 2004]) or 



is to use Eq. (|16p . In this approach, a key parameter is the 
magne tic field, which is likely amplified around t he sho ck 
front (IVink fe Laming) 120031; IBamba et al.ll2003l l2005allbl ; 
lYamazaki et al.l 120041 : lUchivama et all 120071) and may de- 
pend on various physical quantities such as the shock veloc- 
ity, the ambient density, and so on. At present, the evolution 
of the magn etic field is not well understood despite many 
work s (e.g.. iNiemiec et al] 20081: iRiquelme fe Spitkovskvl 
120091 : lOhira et all l2009at iLuo fc Melrose] I2009D . In addi 
tion, the evolution of another parameter rj g is also unknown. 
These facts prevents us from predicting p m rigorously. 

Here we adopt a different phcnomcnological approach 
based on the assumption that you ng SNRs are respo nsible 
for observed CRs below the knee (|Gabici et al.ll2009h . The 
maximum energy cp m is expected to increase up to the knee 
energy (10 15 ' 5 eV) until the end of the free expansion phase, 
isedov, and decreases from that epoch. As seen in § 14.11 p m 

is limited by the escape at t > iscdov, that is p m — pm SC '. 
Then, to reproduce the observed CR spectrum from ~ GeV 
to the knee, we may assume a functional form of 



t 



^Scdov 



— 6.5/a 



(32) 



expected theoretically 



Reville et all 120071: iGiacalone 

2009; lOhira et al.l l2009bl lcT) . This statement is recast if we 
involve recent results of X-ray observations. The precise X- 
ray spectrum of RX J1713.7— 3946 is revealed, which gives 
u sh = 3.3 x 10 8 ?7g /2 cm s" 1 dTanaka et al.ll2008f) Then, 
Eq. (|2"9"]) can be rewritten as (|Yamazaki et al.ll2009f ) 



where cp^ ncc = 10 15 5 eV, so that cpm SC ' 1 = 1 GeV at t = 
10 a ts c dov F° r later convenience, we change variables from 
t to the SNR radius i? s h in order to take x — Rsh, where 
X is a variable appeared in § 13.31 We further assume the 
Jokipil [2003; llnoue et all dynamics of R ah as 



Bel 12 004: 



( 1 



cdov 



(33) 



cp 



(age) 



5 x 10 J S mG 



/ ^agc \ 



where i?s dov is the shock radius at t — isedov Then, we 
obtain 



TeV 



(30) 



^ SC) (^h)=Pkn 



( Rsh 



Hence, in order to obtain cpm < 100 TeV, we need B ^ 
20 fiG in the context of the hadronic scenario of TeV 7- 
rays, which contradicts current estimates of the magnetic 
field. 

A possible solution is to consider the escape-limited 
acceleration. One can find from Eq. (fT7| that if we take 
^ 0.04 pc, the maximum energy becomes 



V -Rscdo 



-6.5/aft 



(34) 



imn{cp^\cp^ c \ 
100 TeV , 



(31) 



which is consistent with the observed gamma-ray spectrum. 
In the following, we consider the model of escape-limited 
acceleration under simple assumptions, estimating the evo- 
lution of the number density and the maximum momentum 
of accelerated particles so as to discuss the spectral index, 

^csc , 

of N csc (p). 
4,2. Evolution ofp m 

Time evolution of the maximum momentum of accelerated 
parti cles, p m , has been so far dis cussed in many contexts 
(e.g., lPtuskin fc Zirakashvilil2003D . One way to estimate p m 



so that we have a = 6.5 jab (see the last paragraph of ^ I3.3|) . 
If we adopt a « 2.5 ad hoc, in addition to b w 2/5 that 
is expected in the Sedov phase, then we find a « 6.5 as 
a phenomenologically required value in the escape-limited 
model. Note that, even if the SNR dynamics is modified by 
the CR escape, b is almost the same as one in the adiabatic 
case and its effect on a is expected to be so small that our 
conclusion is not affected qualitatively. 

4.3. Dynamics of SNR shock waves 

In this subsection, we consider the dynamics of the SNR 
shock in order to estimate the evolution of the normal- 
ization factor of the spectrum K(R s \ l ). A simple treat- 
ment of the dynamics of the SNR shock from the free 
expansion to the adiabatic exp ansion (Sedov) phase has 
been given by several authors (lOstriker fc McKed 119881 : 
iDrurv et all 119891 : iBisnovatvi-Kogan fc Silichl Il995t) . Here 
we extend their method taking account of the cooling by 
CR escape. The total mass of the SNR shock shell is calcu- 
lated as 



M{R ah ) =M ej +47r 



drr 2 p am 



(r) 



(35) 
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where M e j and p am (r) are the ejecta mass and the den- 
sity of ambient gas, respectivel y. The equation of mo- 
tion of the thin shell is giv en by (|Ostriker fc McKedll988t 



tion ot the tmn shell is given by 
iBisnovatvi-Kogan k Silichll 1995T ) 

d(Mu) 



dt 



4nRi h (P in - P am ) 



(36) 



where the gas velocity u is related to the shock velocity 
u s h and the adiabatic index 7 a d as u = 2u s h/(7 a d + 1)- 
Quantities p n and P am are the pressures of the post- 
shock gas and the ambient gas, respectively. For strong 
shocks, one can neglect P am . The explosion energy £ = 
£th + Mu 2 /2 + Q esc consists of the internal energy £ t h = 
47ri?g h i-i n /(3(7 a d — 1)), the kinetic energy, and the energy 
Qcsc which is carried by the escaping CRs. Then Eq. 
is rewritten as 



3 
2.5 
2 

1.5 • 
1 • 

0.5 • 
• 



-0.5 



■ ■ — — 1 

• 
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Without CR cooling 

CR Coolina . 


■ 








■ 













































10"' 



10"' 



10" 

Time[t/t Sed0 „] 



10 l 



10* 



d(Afu) , 3( 7ad - 1) = 



3(7ad - 1) 



dt 



2R 



sli 



R 



sli 



Fig. 1. The index 6 — d In E/d In P s h as a function of time 
for the case PS. The solid and dashed lines are for the 
{£ — Qosc(-Rsh)} .(37) adiabatic expansion and the expansion with CR cooling, 
respectively. See the text for details. 



Since dt = 2dP sh / (7 ad 
becomes 



l)u, the left-hand side of Eq. ([37]) 



7ad 



1 



[Rt h {Muf 



AM stl di? sh 
where u = 6( 7a d — l)/(7ad + !)• Hence we obtain 



(38) 



Wsh(-Rsh) 



7ad + 1 



H 8 h 



2lo 



AP(R sh )R 



sli 



x / dr(8 - Q^y^Mir) 
Jo 



1/2 



(39) 



Once Q esc (r) and p am (r) are given, we can integrate 
Eqs. ([33]) and (]3"9"]> to derive the dynamics of the SNR, 
u s h(-Rsh)- We mainly consider the evolution of w s h in the 
Sedov phase. Hence, in the following, we simply assume 
that /o am is constant with r because even in the case of the 
core-collapse supernova, the wind region is about a few pc 
and has been already passed by the SNR shock until the be- 
ginning of the Sedov phase. The energy carried by escaping 
particles Q C sc is written as 

/•Pknco 

<3csc(iU) = 0(i? sh - i? Sc dov) / c P N csc (p)dp , (40) 

where 9(x) is the Heaviside step function, i.e., Q(x) = 1 
for x > while Q(x) = for x < 0. Substituting Eqs. (g7J , 
and ([M]) into Eq. 1(30]), we obtain 



(s + 0/a ± 2) 

(41) 

(s + /?/« = 2) , 



Qesc(-Rsh) — 7?56(i? s h — -Rsedov) 

l-(fl 3h /flSodov)-° (2 - 3) + 3 

l-(m p c/p kncc ) 2 — P/° 



a 



log(fl sh /-R. S cdov) 

log(p knoc /m p c) 



where 77 is the ratio of the total energy of escaping CRs 
to the explosion energy. Note that 7a d is assumed to be 
constant. Integrating Eq. ([39"]) with the aid of Eqs. |35]) 
and (|4"T]) . the dynamics of the SNR shock is determined. To 
make more realistic discussion, we should perform fluid sim- 
ulations with CR back reaction which is beyond the scope 
of this paper (e.g.. lDrurv et al]ll989t IVolk et al.ll2008f ). 



■ ■ — — 
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CR cooling - 
CR Cooling ■ 





































































10"" 



10"' 



10" 



10 1 
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Fig. 2. The same as Fig. [TJ but for the case PH. 
4.4. Evolution of K 

In this subsection, we discuss the evolution of the normal- 
ization factor, K(R s h), of the spectrum of accelerated parti- 
cles. At present, the injection process for CR acceleration at 
the shock is not well understood. Hence, we consider two 
representative scenarios of the injection process to model 
the amount of the accelerated particles. 

At first, we consider the sam e injection model as that 
of iPtuskin fc ZirakashvilH (|2005l ). The model requires that 
the CR pressure at the shock is proportional to the fluid 



ram pressure, that is, Pgr Pa 
the shock Pgr is given by 



4-7T 



Pm(fl 3 h) 



dpvp 3 f (p) , 



g h . The CR pressure at 



(42) 



where we neglect the contribution of non-relativistic par- 
ticles. Then one can find that Pgr is related to fo(rn p c) 
as 



. /o(m p c) (s>2,PS) 
PCR " { fo(m P c)R:t 2) (*<2,PH) 



(43) 



where we have used the fact that the distribution function 
of CRs at the shock front is essentially /o(p) oc p~ s ~ 2 for 
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p > m p c. Note that Pgr is approximately proportional to 
/o(m p c) when s > 2, because in this case Pgr only weakly 
depends on p m . Hereafter, the cases of s > 2 and s < 2 are 
called PS and PH, respectively. Since Pgr oc p am v% h , the 
normalization factor of the CR spectrum inside the SNR 
shock K(R s h) is calculated as 



K {Rsb) oc P^ h / (m p c) oc 



r E(R sh ) 



i A' 



«(2- 



(PS) 
(PH) 



(44) 



where the mechanical energy of the ejecta is (up to the 
numerical coefficient) 



In summary, the evolution of the normalization factor of 
the spectrum of accelerated particles is given by K (R s h) oc 

rL ^th 

r s (ps) 

(3=1 5 + a{2-s) (PH) . (48) 

I 3+(s low -l)(<5-3)/2 (TL) 

4.5. The spectrum of escaping particles 

We obtain from Eqs. p?|). (|28]) . and (48]), the index of the 
momentum spectrum of escaping particles as 



E(R sh ) 



^PamUsh^sh 



(45) 



If the explosion is adiabatic, E(R s h) is constant with i? s h 

3/2 

during the Sedov phase because u s h oc R sh . However, if 
the modification of SNR dynamics via CR cooling is taken 
into account (§ 14.3ft , then u s h is no longer proportional to 
R~^ 2 , and it is obtained by solving Eq. (|39p . In Figs. [1] 
and [21 we show S = d In Ej d In i? s h as a function of time. 
Figure [1] is for the case PS in which we adopt r\ = 0.2, 
s = 2.23, 7 a d = 5/3, and (3 = -0.2, while Fig. EH is for 
the case PH with parameters 77 = 0.2, 7^ = 4/3, and (3 — 
ot(2 — s) — 0.2. One can find that after t ~ iscdov, <5 is around 
—0.2, so that E{R s \ l ) is approximately given by E{R s \ l ) oc 
i?f h with constant <5 w — O.fl For the case PS, the effect 
of CR cooling appears at the late time because low-energy 
CRs have a large fraction of the total energy of the escaping 
CRs. On the other hand, for the case PH, the effect emerges 
at an earlier epoch. 

Next, we consider the thermal leakage (TL) model 
(jMalkov fc Volkl fl995). This model requires the continu- 
ity of the distribution function fo(p) to the downstream 
Maxwelian at the injection momentum pj n j, namely 



/oOinj) OC Pam/pfnj ■ 

Using this fact, we derive 
K(R sh ) oc R^ h f (m p c) oc p am R^ h p 



(46) 



(47) 



where si ow is the spectral index in the non-relativistic 
regime, that is, fo(p) oc p~ s i°w-2 f or ^.^ < p < mpCi \y e 
further assume the constant p am and p lni oc u s h- Then we 
obtain K(R sh ) oc i? 3+( s >°w-i)('5-3)/2^ go that ^ = 3+( < Siow __ 
1)(5 — 3)/2. Although the value of s\ ow is uncertain, one can 
expect s\ ow > 2. In the test particle approximation, one 
derives si ow = s > 2. If the nonlinear model of shock accel- 
eration is considered, the spe ctrum fp(p) is softer than p 
in the non-relativistic regime (jBerezhko &: Ellisonlll999f) . In 
particular, if s\ ow w 3, then f3 ~ (5, so that the value of <5 
becomes similar to those of PS case. 



1 We search for approximate solutions to Eqs. (|39J1 . (|4ip and 
which determine S (and /3) for given parameters a, s, rj and 
7 a d. The procedure is as follows. For a trial value of /3 which is 
assumed to be constant, we solve Eqs. (|39[) an d (|41[) to obtain 
u s h (-Rsh) so as to calculate E and <5 as functions of i? s h- It is 
found that 5 is almost time-independent after isedov, hence we 
can derive its average value in the epoch t > ts c dov Then, we 
update the value of /? using Eq. (|48[) with the average value of 
8. We repeat this procedure until the iteration converges. 



C s + S/a 
= \ 2 + 5/a 
{ s + {3 + (sio 



(PS) 
(PH) 

l)(*-3)/2}/a (TL) 



(49) 



In the following of this section, we discuss which injection 
model is suitable to make the galactic CR spectrum ob- 
served at Earth. We adopt a — 6.5 as a typical value (see 
ij 14. 2p and s esc = 2.2 to make the Galactic CR spectrum. In 
this paper, we assume s csc « 2.2-2.4, where we believe that 
the deviation from 2.0 is significant. This would be tested 
not only from Galactic CR observations but also observa- 
tions of extragalactic galaxies. Note that, even if we adopt 
s C sc = 2.4, our results are qualitatively unchanged. 

In the case PS, s csc is smaller than s because S 0, 
so that the model predicts the harder spectrum of es- 
caping particles than that of the source. However, since 
5/a « 0.03, the difference is small. In order to reproduce 
the observed Galactic CR spectrum s osc , the source spec- 
trum should be s ~ s csc « 2.2. Hence, the PS model re- 
quires s > 2 at the source. This condition is satisfied if 
we consider th e diffusive shock acceleration at moderate 
Mach number (|Fuiita et al. 2009) . It is also possible to de- 
rive s > 2 if we c onsider the effects of neutral particles 
(|Ohira et al.ll2009bl lcM. 

In the case PH, because the value of S/a is small, s esc 
always near 2, which is the value predicted by the diffusive 
shock acceleration theory in the strong shock, test-particle 
limit. In particular, if S = — indeed, even in the test- 
particle limit where the cooling via CR escape can be ne- 
glected, the value of 5 is not exactly zero unless t ^ ^Sedov 
(see Fig. [H — then, o ne ca n find s esc = 2. This is what 
iBerezhko fc Krvmskvl (<1998f ) and iPtuskin fc ZirakashvilH 
( 20051 ) showed. Note that from Fig. El 5 is negative for a 
long time, so that s csc < 2. Therefore, this model PH can- 
not reproduce the observed Galactic CR spectrum at the 
Earth. 

In the case TL, neglecting <5, we obtain 



s csc ~ S ~t~ ^low) 



(50) 



Hence if s\ ow < 3, then s esc > s. For the test particle ac- 
celeration, we have s = si ow , so that the source spectrum 
should be s « 2.2 in order to obtain s osc = 2.4. For the 
nonlinear acceleration, one typically expects s\ ow > 2 and 
s < 2, so that s csc < 29/13 « 2.23. Therefore, the non- 
linear model seems to explain the inferred source spectrum 
marginally. Possible rigorous determination of the value of 
s csc may give a constraint on the theory of nonlinear accel- 
eration. 
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5. Application to AGN Cocoon Shocks 

In this section, taking account of a constraint derived 
from the spectrum at the Earth, we study the ori- 
gin of CRs with energies larger than ~ 10 175 eV and 
the acceleration mechanism of them at AGNs. There 
are many wo rks which discuss UHECR pro d uction in 
AGNs ( e.g., iBiermann fc Strittmatterl 119871 : iTakaharal 



1990; Rachcn fc Biermannl Il993l; iBerezinskv et al.l [2006: 



Berezhkol 120081: iPe'er et al.ll2009D . Many of them focus on 
UHECR acceleration in radio galaxies including Fanaroff- 
Riley (FR) I and II galaxies, which typically have powerful 
jets. In the context of DSA, one can basically suppose three 
acceleration zones; internal shocks in jets, hot spots, and 
cocoon shocks. The former two are most widely discussed 
scenarios but the detailed study of DSA at such mildly rel- 
ativistic shocks has not yet been achieved. In this section, 
we concentrate on the cocoon shock scenario proposed by 
Bcrezhko (2008), where the non-relativistic DSA theory can 
be applied. 

In this scenario, extragalactic CRs with energies larger 
than the second knee (~ 10 175 eV) may be accelerated 
at the outer cocoon shock running into the intergalactic 
medium (IGM). As powerful jets penetrating into a uni- 
form ambient medium with a density Pa, m , the heads of the 
cocoon advances into the IGM with a velocity u^- At the 
same time, the cocoon expands sideways with a velocity u c . 
Since the typical cocoon shock is non-relativistic, we apply 
the escape-limited model considered in previous sections. 
Although we hereafter focus on this scenario, note that it 
is very uncertain whether the efficient acceleration occurs 
there since observed non-thermal emission is much weaker 
than that from hot spots and lobes. 

In the following, we investigate whether the CR spec- 
trum above the second knee can be explained by the AGN 
cocoon shock scenario with the same parameters for young 
SNRs explaining the CR spectrum below the knee, which 
have been discussed in § 2J Similar to the previous cal- 
culations in § ID we hereafter calculate the values of a = 

d In pm /d In t and j3 = d In K/d In t in order to derive the 
spectral index, s esc . Here we adopt x = t, where x is a 
variable appeared in § 13.31 

First, let us consider evolution of the maximum mo- 
mentum, pm C ' (i) , in a phenomenological way. In the young 
SNR case (§ 14. 2|) , we have phenomenologically expected 
p m oc t~~ . Then, by using the Sedov- Taylor solution 
(u s h oc £~ 3 / 5 and R s h oc t 2 / 5 ), we can easily obtain 



P 



(° sc )(i)(xu sh 3 JE 



■sh , (51) 

where c is a phenomenologically introduced parameter since 

one can expect I oc i?° h in general. Note that Pm oc B£u B h, 
however, B does not depend on R s h but on it s h as long as B 
is generated by plasma instabilities such as CR streaming 
instabilities. In the cocoon shock scenario, u s h and i? s h are 
replaced with u c and i? c , respectively. 

In order to obtain the value of a — dhipm^ /dlni, the 
dynamics of the AGN cocoon is necessary. A simple con- 
sideration of the cocoon dynamics for the constant density 
IGM tells us that is almost time-independen t and u c 
evolves as u c oc t^ 1 / 2 (Be gelman fc Cioffil[T989l) . so that 
the cocoon radius evolves as R c oc u c t oc t 1 / 2 and the jet 
radius evolves as i?h oc Uht oc t. Then, we obtain a = i 3 - 
and a — 2 for c — and c = 1, respectively. 



Next, let us consider the time dependence of the normal- 
ization factor of the spectrum of accelerated CRs, K(t) oc 

. The volume of the acceleration region swept by the co- 
coon shock is ~ Au c t, where A is the total area of the 
shock surface. If we assume the elliptical shape of the co- 
coor0, then A oc RcU^t oc i 3 / 2 , so that Au c t oc t 2 . In the 
cases PS and PH, /o(m p c) is related to Pgr- The depen- 
dence of Pgr is written as Pgr oc Pam^c x '■> where we 
neglect for simplicity the evolut ion of the accele ration effi- 
ciency which was considered in iBerez hkol (I2008h . Then, 111 
the case PS, fo(m p c) oc Pgr leads to 

K(t) oc (Au c t) x / (m p c) oc i 1 , (52) 
while in the case PH, /o(m p c) oc Pcr^ 2-3 -' leads to 
K(t) oc (Au c t) x / (m p c) oc t 1+3t( - 2 ~ s) . (53) 

Finally, in the case of model TL (see § 14. 4p . /o(m p c) oc 
PamP^r" 1 leads to 



K(t) oc (Au c t) x /o(TOpc) oc t( 5 - s '=w)/2 



(54) 



where we assume the constant IGM densitjfj Therefore, we 
obtain 



. r 1 (ps) 

(3= { 1 + 5(2 - s) (PH) 
I (5 - s low )/2 (TL) 



(55) 



By using the above results, we can obtain the spectrum 
of escaping particles. Here, we assume a — 2 expected for 
c = 1 (I oc R c ). Then from Eq. (|2"5j). we derive 



f s + 1/2 (PS) 
= 5/2 (PH) 
{ s + (5 - s low )/4 (TL) 



(56) 



Interestingly, all the three cases (PS, PH and TL) lead 
to the source spectral index s csc « 2.4-2.7 which is re- 
quired in the proton dip model and in the other extra- 
galactic scenarios. In the case of model PS, s = 2.2 leads to 
Scsc = s + 1/2 = 2.7. In the case of model PH, s < 2 leads 
to s csc = 2+1/2 = 2.5. Note that in this case, the value 
of s esc does not depend on s but on a generally. Finally in 
the case of model TL, s = s\ ow = 2 (based on the test par- 
ticle acceleration) leads to s csc = 2 + 3/4 = 2.75, which is 
also consistent with the observed UHECR spectrum in the 
proton dip model unless the redshift evolution of UHECR 
sources are fast. If the effect of nonlinear acceleration is 
prominent, s < 2 and s\ ow > 2 can give harder indices of 

°csc 

< 2.75. 



2 iBerezhkol (|2008T ) assumed a kind of spherical cocoon, i.e., 
A oc R c u c t, and used K(t) oc p am E|! which is similar to the 
relation obtained in the case of model PS. However, when the 
cocoon becomes spherical, u c and evolve according to the 
adiabatic solution as in the Sedov- Taylor solution for SNRs 
(|Begelman fc Ciom|[l989l : iFuiita et al.ll2007h . 

3 However, it may not be the case since the IGM density may 
drop with the distance from the nucleus. Then we need to use 
different adiabatic solutions for u c and R c . 
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6. Summary and Discussions 

In this paper, we have investigated the escape-limited 
model of CR acceleration, in which the maximum energy 
of CRs of an accelerator is limited by the escape from the 
acceleration site. The typical energy of escaping CRs de- 
creases as the shock decelerates because the diffusion length 
becomes longer. After revisiting the escape-limited model 
and reconsider its detail more generally, we have derived 
a simple relation between the spectrum of escaping parti- 
cles and one in the accelerator. Then, using the obtained 
relation, we have discussed which model of injection is po- 
tentially suitable to make the Galactic and extragalactic 
CRs observed at the Earth. 

For young SNRs, we have considered the shock evolution 
with cooling by escaping CRs and those spectra for the 
three injection models. As a result, we have found that in 
the case PH, it is difficult to satisfy the condition for the 
source spectrum of Galactic CRs (s osc « 2.2-2.4). On the 
other hand, s csc 2.2-2.4 can be achieved in cases of PS 
and TL. 

We have also applied our escape-limited model to AGN 
cocoon shocks as well as young SNRs. This model is just 
one of the various candidates proposed so far, even if AGNs 
are UHECR accelerators. Nevertheless, it is interesting that 
the young SNR and the AGN cocoon shock scenarios can 
explain the Galactic and extragalactic cosmic rays observed 
at the Earth in the same picture for all the three injection 
models if we accept the proton-dip model inferring s csc ~ 
2.4-2.7. Whether the proton dip model is real or not can 
also be tested by future UHECR and high-energy neutrino 
observations. In this paper, we have focused on the proton 
case. Obviously, heavier nuclei become important above the 
knee so that we need to take into account them in order to 
explain the CR spectrum over the whole energy range. We 
can also apply the escape-limited model to heavy nuclei 
CRs for this purpose, although it is beyond the scope of 
this paper. 

We point out a potential problem for the magnetic field 
amplification in the escape-limited model. In the case of 
young SNRs, we have determined the evolution of the max- 
imum energy in the phcnomcnological way, and adopted 

10 + 2e 

a 6.5. Using Eqs. (fTrj)) and (|5"Tj) . we obtain B cx u sh 3 , 
where i cx R^ h . The same result is obtained for AGN co- 
coon shocks because we have considered the case in which 
the same parameters describe both the young SNR shocks 
and the AGN cocoon shocks. In particular, for c = 1 as 



is used in Eq. (31), we obtain B 



oc u" 



which means 



that B rapidly decreases with radius (or time). In prin- 
ciple, both the dependence of B on u s h and the value of c 
can be determined theoretically, and then the evolution of 
the maximum energy should be predicted. Some previous 
works are based on theoretical arguments on the magnetic 
field evolution (e.g., 
l2008tlCaprioli et al 



Bell2004: Pcllt ier et al.ll2?jol [Berczhko 



2009), which seem to be different from 



our phenomenological one. At present, the mechanisms of 
the particle acceleration and the magnetic field amplifica- 
tion are still highl y uncertain despite of many theoretical 
efforts ( e.g., Niemiec et all 20081: Ricmelme fc Spitkovskvl 



120091: lOhira et al.ll2009aULuo k Melrose! l2009f ). Hence, we 
expect that further theoretical and observational studies 
might reveal this discrepancy or exclude the possibility of 
escape-limited acceleration in the future. 



In this paper, we have mainly considered spectra of dis- 
persed CRs around young SNRs and AGN cocoon shocks. 
However, applications to other astrophysical objects are, 
of course, possible. For example, the old SNRs detected 
by Fermi LAT, s uch as W28, W44, W51 and IC443 
(jAbdo et al.l l2009allH ) have been of great interest because 
they likely generate escaping CRs. In fact, the number of 
CRs around such old SNRs is likely to decrease with time 
or the shock radius, that is (3 < while a > 0. For ex- 
ample, when we consider the dy namics of an old SNR , we 
have u s h oc i~ 2 / 3 cx R~ h 2 (e.g.. lYamazaki et al.ll2006T ). so 
that we have E cx u 2 R^ cx R~^, i.e., 5 = — 1 < 0. On 
the other hand, the value of a may be different from 6.5, 
which could be attributed to various complications such 
as the interaction with the dense molecular cloud, and so 
on. For example, for the maximum hardening case, that 
is, s csc = s + (s — l)(cr _1 — l)/a (see Appendix A. 2), we 
find Scsc ~ 1-5 when a w 1 and s ss 2.5 where we assume 
s = (cr + 2) /(a - 1). This might b e the case for old SNRs 
such as W51C (jAbdo et al.ll2009ri l. In addition, the maxi- 
mum energy may be rather small for the old SNRs, so that 
the spectrum above pm°' would be suppressed. The spec- 
trum of high-energy gamma rays might give us important 
information on both the acceleration and escape processes 
of CRs with energies much lower than the knee energy. 
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Appendix A: Effect of Adiabatic Loss on the 
Spectrum of Escaping CRs 

The instantaneous spectrum of CRs escaping from the ac- 
celeration site when the shock radius R s h is obtain ed from 
the equation (23) of iPtuskin k. ZirakashvilH (|2005ft . 



d/e; 



dR 



sli 



1 d/ csc 
u s h dt 
AirS (p — p 



(l-L-Ss. 

\ x a 2 



shl 



Rsh ( dp 

drr 2 f(p m ,r) 



(A.l) 



where a and £ cr are the total shock compression ratio and 
the ratio of the CR pressure to the shock ram pressure. The 
spectrum and the maximum momentum at the shock are 
assumed as 



f(p,R sh ) = Ap- s - 2 R^ 3 , 

Pm(Rsh) = Rs\? , 



(A.2) 
(A.3) 



where A is a normalization factor of the distribution 
function, and p and p m are normalized by pknco while 
i? s h is normalized by i?scdov The distribution function of 
CRs at r < f? s h can be fo und by the same method as 
IPtuskin fc Zirakashvilil(|2005l ). The fluid velocity at r < R sh 
is 



tt(r '* )= ( 1 -£)idk ,, * w ' 

and we solve the following equation 
df ,. M df 1 d . , , Of 

Then one can get the following solution. 

= A n s-2 r a( S +l3-l)-s-2 R - (s+(i- 1) (a- 1) 



f(p,r)=Ap- s -'r 



sh 



(A.4) 



(A.5) 



(A.6) 



Therefore the total spectrum of escaping CRs is 



fesc 



d/e: 



-dR sh 



di? sh 

— fesc, surface H~ fese, inside ; 



(A.7) 



where / esc .surface IS the firs t term Of Ecj. (|A. 1[) and /esc, inside 

is the second term of Eq. (| A. 1|) . These are given by 



At: A 1 



l _ C„ 



csc, surface 



_2_ -s-2-/3/a 

l J J 



(A. 



/ csc,insi. 



dc 



= At: A [a - 1 + - 



x / dRSip-p^mp-i-sR-^-^-v- 1 

Jo 

r R 

x / drr a( - s+fJ - 1 ^ s . (A.9) 
Jo 



The dr integral of Eq. (|A.9[) is approximately given by 

( j rr ,<r(s+/3-l)-s 

j ^( s+/3 -t)- s+1 [fJ >{s _ 1){ L_ l)] 

where Ro is the radius at which /3f roc (> (s — 1) — l)) 
changes to (3 (< (s — 1) — l)). Here, we assume that 
f3 changed when i? s h = 1, that is, t = t S edov Then from 
Eq. jXl, R =Rfc< r ~ 1 . 
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AJ./3>(s-l)(i-l) 

In this case, the dr integral of Eq. (|A.9[) is dominated by 
the outer region, that is, spectrum of escaping particles 
does not depend on the past acceleration history. / eS c, inside 
is calculated as 

47rA(a-l + I) 

J esc. inside 



1 - s + ct(s + j3 - 1) 

x /** dRS( P - Pm (i?))p m 1 - S ^- 1 
Jo 

= ^ ^4±i -p— (A.ii) 

a l-s|d(sl/1-l) 

Hence, as long as /? > (s — 1) (a -1 — l), the energy spec- 
trum index of escaping CRs is s csc = s + (3 /a. 



A.2. /3<( a _i)(I-i) 

In this case, the dr integral of Eq. (|A.9[) is dominated by 
the inner region, that is, the spectrum of escaping parti- 
cles depends on the past acceleration history. Especially, 
only the acceleration at t = i S cdov is important, and then 
/esc,insidc is calculated as 



/esc, inside 



4irA(a - 1 



1 - s + cr(s + /3 froc - 1) 



X 



/ " dRS{p - Pm (i?)) Pm 1 - s i?( s - 1 )^- 1 )- 1 

JRo 



= 4ttA + - p - s - 2 - (s _ 1)( X- 1)/a ^ 

a 1 - s + cr(s + /3f roo — 1) 

(A.12) 

Because / C sc,insidc is softer than / eS c, SU rface (at p < p m ), 
/escansidc is larger than f esc<SUI { ace . Hence the spectral index 
of escaping CRs is s csc = s+ (s — l)(o" _1 — 1)/q and does 
not depend on 0. The largest possible hardening from the 
spectral index of the acceleration site As = s — s csc is 

s- 1 / 1 

As = 1 - - 

a \ er 

Qf(s + 2J 



where we assume the relation between the spectrum index 
and the compression ratio, s = (cr + 2)/(<r — 1). 



